Static versus dynamic heterogeneities in the D = 3 Edwards- Anderson-Ising spin glass 
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We numerically study the aging properties of the dynamical heterogeneities in the Ising spin 
glass. We find that a phase transition takes place during the aging process. Statics-dynamics 
correspondence implies that systems of finite size in equilibrium have static heterogeneities that 
obey Finite-Size Scaling, thus signaling an analogous phase transition in the thermodynamical limit. 
We compute the critical exponents and the transition point in the equilibrium setting, and use them 
to show that aging in dynamic heterogeneities can be described by a Finite-Time Scaling Ansatz, 
with potential implications for experimental work. 
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Spin glasses, fragile molecular glasses, polymers, col- 
loids, and many other materials display a dramatic in- 
crease of characteristic times when cooled down to their 
glass temperature, Tg This is probably due to the col- 
lective movements of an increasing number of elements 
in the system, with a (free) energy barrier growing with 
the size of the cooperative regions Q (the cooperative 
regions become larger as the temperature gets closer to 
Tg). Experimentally, one can get the fingerprints of these 
movements by observing dynamical heterogeneities or 
non-linear susceptibilities [4L 

Below Tg, aging appears 5[- Consider a rapid quench 
from a high temperature to the working temperature T 
(T < Tg), where the system is left to equilibrate for time 
tw and probed at a later time t + t^. One finds that the 
response functions (e.g., magnetic susceptibility) depend 
on t/ti^, with /i « 1 j5|-|7|. The age of the glass, t^, 
remains the relevant time scale even for ~ days. 

Dynamical heterogeneities age as well, as found numer- 
ically in their characteristic length ((t,tvt) [1, 0]. Recent 
measurements of aging correlation and response functions 
with space-time resolution suggest that C(t,iw) will 
soon be experimentally investigated. Characterizing ag- 
ing for ({t,tv,) is our main concern here. 

We focus on spin glasses, an easier case for a number 
of reasons: (i) the sluggish dynamics is known to be due 
to a thermodynamic phase transition at Tc = Tg [ll- 13|: 
(ii) the size of the glassy magnetic domains, ^(iw), is ex- 
perimentally accessible [3, [T^ 100 lattice spacings 
at T ~ Tc [l^ . larger than comparable measurements 



for structural glasses Q); (hi) ^(t^) oc tl/^'''^\ z{T) « 
6.9Tc/T suggests that free-energy barriers grow in spin 
glasses as ~ log^(tw), rather than with a power law as 
in fragile glasses; (iv) equilibrium physics is known to 
rule nonequilibrium dynamics [16|. A quantitative cor- 
respondence exists between equilibrium and nonequilib- 
rium spatial correlation functions 17, [l3| (equilibrium 



on systems of size L matches nonequilibrium at time tvj). 
Finally, the Janus dedicated computer [l^ allows us to 
simulate nonequilibrium dynamics from picoseconds to 
a tenth of a second [^, [l3l , and to compute equilibrium 
correlation functions on lattices as large as T = 32, and 
temperatures as low as 0.64 Tc fisj . 

In this paper we show that a phase transition occurs 
in the aging dynamic heterogeneities. As time t pro- 
ceeds, when the spin correlation function C{t,t^) (see 
below) becomes smaller than the spin glass order param- 
eter qea, the length scale of the dynamic heterogeneities 
C(<,iw) diverges in the limit of large t^. We use the 
statics-dynamics correspondence to investigate this phase 
transition in the equilibrium setting, focusing on spatial 
correlation functions (siatic heterogeneities). Finite-Size 
Scaling (ESS) yields an accurate estimate of (;ea (some- 
thing never achieved before for a spin glass) as well as 
of the relevant critical exponents. Back to nonequilib- 
rium, aging turns out to be amazingly well described by 
a Finite-Time Scaling Ansatz, with critical parameters 
taken verbatim from the equilibrium computation. 

We consider the Edwards- Anderson model on a cu- 
bic lattice of size L (volume V = Li^), with periodic 



2 



boundary conditions, at T = 0.64Tc. We use Ising spins, 
Sx = ±1, and binary nearest-neighbor couplings. The av- 
erage over the quenched disorder, denoted by an overhne, 
is taken after the thermal average (...). We consider two 
clones of the system, {sx\sx^} evolving independently 
under the same set of coupling constants, and taken at 
the same time t^. The replica field is = s^^s^x^ and 
the spin overlap is its spatial average q = qx/V- See 
Rcfs. and for full details of our equilibrium and 
nonequilibrium simulations. 

Out of equilibrium, correlation functions depend either 
on a single time ty^,, or on t and t^. Let Cx{t,t^) = 
Sx(t + tw)sx(iw)- The spin correlator, see Fig.[T]-top, is 

^w) = ^ {Cx{t,t^)) , Cit) = Cit, U=(X^). (1) 

X 

Naive aging is approximatively valid: for finite t^, 
C{t, tv,) decays for long t, but the decay slows down with 
increasing In fact, there is an enveloping curve C{t) 
with a non-zero limiting value, the order parameter qea- 
The lack of a reliable parameterization of C{t) precludes 
a controlled extrapolation of (/ea, in contrast with the 
equilibrium computation shown below. 

As for space dependencies, we consider C4(r,iw) = 
Yl^{9=oitw)qx+rit^))/V. Using integral estimators d, 
17[ we extract the coherence length ^(iw), the size of 
regions where the two clones of the system are simi- 
lar. Yet, to learn about heterogeneities on the dynam- 
ics probed at time t + tvr, at distance r, we consider 

C22(r,i,tw) = Ex i'^dt, U)Cx+r(.t, iw)) " C^it, U)/V . 

Using an integral estimator [9|, we extract from C22 the 
correlation- length C,{t,t^), the characteristic length for 
heterogeneities, displayed in the central panel of Fig. [T] 
We replace t with C{t,tw) (l^l, as independent variable. 
For large C, C(C':*w) reaches a ^w-independent value, 
which increases when C decreases. On the other hand, 
for small C, C(C, iw) grows strongly with t^j- Clearly, 
something happens when C goes through some special 
value f^EA and we intend to exploit the statics-dynamics 
correspondence to clarify it. 

How does all this appear from an equilibrium view- 
point? In the limit of large system size L, the probabil- 
ity density function for q, P[q) = Pj{q), has two Dirac's 
delta contributions of equal weight at g = i^EA- Replica 
Symmetry Breaking (RSB) theory [2l| predicts that P{q) 
has a support for \q\ < qv A , w hile droplet theory expects 
no support in that region [22|. 

Our approach focuses on the study of equilibrium con- 
nected correlation functions [23| , regarded as a function 



of the spin overlap q. Varying q at fixed T a phase transi- 
tion is encountered for q = qEA- As in jl8[, our conditional 
correlation function at fixed q = c, C4(r|c), is obtained 
as a quoti ent of the convolutions of (qxqx+rS{q — c)) and 
{6{q — c)) with a Gaussian of width 1/V. This combines 
0{\/V) levels, thus smoothing the comb-like P{q) p3 |. 
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FIG. 1; (color online) Top: C{t,t^), Eq. ((T)), as a function 
of t for = 4\ i = 3, .., 16 (lines, grows from bottom 
to top). We also plot C{t) (points) and our result for qea 
from Eq. (horizontal line). Center: Correlation length 
^(C, fw) as a function of (same values of as in the top 
panel). Bottom: Finite-time coherence length ^(t^) against 
the finite-size coherence length at q = 0, (,{L). The results 
are compatible with ^(tw) = 1.48^(1/) (straight line). 



It has been recently found [17|, [18| that the equilib- 
rium C4{r\q), computed in a system of size L, accurately 
matches the nonequilibrium C22{r,C,t^) if one chooses 
time t such that C{t,t^) ~ q and time such that 
L « 3.7$(tw) (at least at T = 0.64Tc). It is tempting 
to assume that the correspondence will become exact in 
the limit of large L and t^. In Fig. [l]-bottom we show 
an example of this correspondence in the limit C — >■ 0. 

To proceed with the equilibrium analysis, we observe 
that C4{r\q) tends to q^ for large r. In a finite sys- 
tem, one needs to perform a subtraction that compli- 
cates the analysis [23[. We instead consider the Fourier 
transform at wave vector fc, C4{k\q) = Er G^''"^C4{r\q), 
blind to a constant subtraction for fc 7^ 0. Defining 
fcmin = (Stt/L, 0, 0) (or permutations), wc have 

= C4(femin|g). (2) 

For T < Tc and \q\ < qEAi one expects that 
A, 



Ciir\q) ^ q' 



C4ik\q) cx ^"(9'-^ + 



(3) 



(scaling in Fourier space holds only if 9{q) < D). The 
dots in ([3]) stand for scaling corrections, subleading in 



3 



the limit of large r (or small k). On the other hand [2 

1 



Ciik\q^ > q^j^) cx 



fc2 



(4) 

The correlation length diverges when jg] — >■ qea from 
above, oc (g^ — Qea)"''- principle, i> is different 

from the thermal critical exponent at . We note as well 
the scaling law (27| 



[qEA) 



2/v. 



(5) 



The theories for the spin-glass phase differ in the pre- 
cise form of 0{q), but agree that a crossover can be de- 
tected in Fq for finite L. 
\q\ < (jEA, while Fq I 
For large L, ^^=°° cx (g^ - ^Ia)" 



Indeed, Fq ~ L'^-^'^i) for 



> gEA, see Eqs. P4| . 
and the crossover be- 
comes a phase transition. FSS tells us that, see e.g. [25| . 



for 

,2 



-D-e(gEA) 



G(il/'^(g-gEA)), 



(6) 



up to scaling corrections. G is a scaling function. 

We have exploited Eq. ([6]) in the following non- 
standard way. We focus on quantities depending on the 
continuous parameter y (e ~ D — 6{q-EA))- 



Fq/iy ^L'-yG{L'/'{q-qEA)) 



(7) 



When y is smaller than D — 6{0), Fq/L^ vanishes in the 
large-L limit for |g| > ^eA: while it diverges for |g| < q^A- 
Hence, fixing y, the curves for pairs of lattices {L,2L), 
will cross at a point qL,y, see Fig. [21 To leading order in 
L~^, the crossing point approaches qea for large L as 



QEA 



AyL 



-1/0 



G(o) ~ 1 
G'(0) 2i/'> - 



(8) 



Note that the amplitude Ay changes sign at y~e. 

We could use a fit to Eq. ([S]) in order to obtain the 
order parameter, but, for a fixed y, there are only three 
crossing points (L — 8, 12, 16) for three parameters (zero 
degrees of freedom). Fortunately, one can extract more 
information from the data by computing the crossings 
for several y values and performing a joint fit, sharing 
gEA and l/v. Since these additional crossing points are 
not statistically independent, this procedure requires a 
proper consideration of the cross-correlations. This can 
be achieved by computing the fit goodness estimator 
with the full covariance matrix for the qL,y The number 
of y values considered is a compromise between adding 
more degrees of freedom and keeping the covariance ma- 
trix invertible. We have chosen 9 values of y obtaining 
= 18.9, reasonable for a fit with 16 degrees of freedom 
(see Fig.O. The result is 



gEA = 0.52(3), l/;> = 0.39(5). 



(9) 



These numbers are remarkably stable to variations in 
the set of y values. Also, removing the L = 8 data for 




FIG. 2: (color online) Fq/L^, Eq. 0, against for all our L 
values at T = 0.703 ~ 0.64rc, both for the free-field scaling 
(F,(fe) oc i.e., y = 2) and for y = 2.35 ^ D - d{qEA). 

The insets are closeups of the crossing region. 
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FIG. 3: (color online) Crossing points qL,y for Fq/V com- 
puted from pairs of lattices (L, 2L), versus L-i/^, for different 
y values at r = 0.703. The continuous lines are fits to Eq. ((8)1, 
constrained to yield common values for (/ea and l/i>. 
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FIG. 4: (color online) Dimensionless ratio C,{C,t^)/£,{t^ 
against the scaling variable (C^ — q^A)i{^v,Y^'' ■ 
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y = 2.3, 2.4 (the outliers in Fig. |3]) shifts our results by 
one fifth of the error bars. Note as well that the slope 
Ay changes sign at y « 2.35. Hence, O^qea) ~ 0.65, in 
reasonable agreement with Eq. ([S]). 

The value of gEA computed above should be the same 
as the large-L extrapolation of the position of the peak 
in P{q): a fit qEA{L) = Qea + AL^^/^, with v from (O, 
yields ^ea = 0.54(3) [H. 

We are finally ready to discuss aging in the dynamic 
heterogeneities. The statics-dynamics correspondence 
suggests that it will take the form of a Finite- Time Scal- 
ing (FTS) Ansatz, similar to Eq. (jG)), in which £,{t^) plays 
the role of L. Since it is a length, C(C, tw) should have 
the same scaling dimensions of ■CC^w)- Setting short-time 
corrections aside, Fig.|4]shows indeed that ^(C, tw)/^(iw) 
behaves function of (C^ - 9|a)?(^w)^/^- FTS also 
provides a natural explanation for the extremely small 
exponents found in i-extrapolations for C{t) j28| . 

In summary, we have studied aging properties in glassy 
dynamic heterogeneities for the Ising spin glass, charac- 
terized through their characteristic length ({t, tv,)- Aging 
takes the form of a Finite-Time Scaling Ansatz, which 
describes the crossover from a regime where C(t,tw) is 
of order one, to a regime where it is of order ^(iw), the 
coherence length yielding the size of the glassy magnetic 
domains. In the limit of an infinite waiting time, the 
crossover evolves into a phase transition. We have prof- 
ited from the statics-dynamics correspondence [13 
to study this phase transition via equilibrium spatial 
correlation functions, thus obtaining the critical expo- 
nents and, for the first time, the spin-glass order param- 
eter. These critical parameters, taken verbatim, describe 
our noncquilibrium data. For a discussion of the mode- 
coupling transition in glass-forming liquids in a similar 
vein see [26| . 
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If we add an interaction h q to the Hamiltonian, 
(i) the correlation length diverges as ^(/i) oc h~^^; (ii) 
h oc [q^A{h) - gEA(0)]i''[i-''-(°-*('''=A", since dgEA/d/icx 
^D-fl(gEA). and (iii) e{qEA)^2{D~i^-') (see, e.g., [H). 



Observe that ^ + ^ — D. 



with a ~ 0.05. Clearly 

qBA(xat)-'^' [C>) = 
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[28] In [g|, we find C{t) ~ q^A oc T 
enough, FTS implies that C{t) 
(^{t,t^ = 00)]. Now, full-aging (as well as empirical evi 
dence [i|) suggest that ({t) oc fT/[Tez(To)]^ j^^^ ^^^^-^ 
Therefore a = T/[Tci>z{Tc)] ~ 0.04. 



